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Abstract: We develop a systematic approach to describe the finite lifetime effects in 
the threshold production of top quark-antiquark pairs. It is based on the nonrelativis- 
tic effective field theory with an additional scale p l l 2 mt characterizing the dynamics of 
the top-quark decay, which involves a new expansion parameter p = 1 — mw/mt- 0m 
method naturally resolves the problem of spurious divergences in the analysis of the unsta- 
ble top production. Within this framework we compute the next-to-leading nonresonant 
contribution to the total cross section of the top quark-antiquark threshold production in 
electron-positron annihilation through high-order expansion in p and confirm the recently 
obtained result. We extend the analysis to the next-to-next-to- leading 0{a s ) nonreso- 
nant contribution which is derived in the leading order in p. The dominant nonresonant 
contribution to the top-antitop threshold production in hadronic collisions is also obtained. 



Keywords: [Heavy Quark Physics, Standard Model, LEP HERA and SLC Physics| . 



Contents 

[l]. Introduction |l| 

p|. Finite width effect beyond the complex energy shift |2] 

|3|. Relativistic corrections |7| 

Strong coupling corrections 

[H]. Eliminating the spurious divergences |9| 

|6|. Results [11 

0. Measured cross section and invariant mass cuts 13 

|8|. Summary 15 



1. Introduction 



The threshold production of top quark-antiquark pairs at a future linear collider may pro- 
vide us with the most accurate information on the top-quark mass and couplings crucial for 
our understanding of electroweak symmetry breaking and mass generation mechanism fij . 
Due to renormalization group suppression of the strong coupling the nonrelativistic top- 
antitop pair is the cleanest quarkonium-like system. Its theoretical description is entirely 
based on the first principles of QCD and is an ideal laboratory to develop the nonrelativistic 
effective field theory approach. It is not surprising that since the pioneering papers f2), ||, Q 
the top-antitop threshold production remaines in the focus of theoretical research for over 
two decades. A significant progress has been achieved in the analysis of the higher order 
perturbative and relativistic corrections in the strong coupling constant a s and the heavy- 
quark velocity v. Sizable next-to-next-to-leading order (NNLO) corrections to the total 
cross section have been found by several groups ||, ^, [?], ||, [9| that stimulated the study 
of the higher orders of perturbation theory. Currently a bulk of the third order correc- 
tions is available 0, [U], [H| |l|, [U], [U| [17|, [l|, |l|, |0|, [21], §2], |24j with only a few 
Wilson coefficients still missing, and the N 3 LO analysis is likely to be completed in the 
foreseeable future. Moreover the higher order logarithmically enhanced corrections have 



been resummed through the effective theory renormalization group methods [25, 26, 27]. 

At the same time much less attention has been paid to the analysis of the effects related 
to the instability of the top quark pq , |29| , [30| , 31 1 . The width of the electroweak top-quark 
decay t — > W + b, Tt ~ 1.5 GeV, is comparable to the binding energy of a would-be toponium 
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ground state and has a dramatic effect on the threshold production. It serves as an infrared 
cutoff, which makes the process perturbative in the whole threshold region, and smears out 
all the Coulomb-like resonances below the threshold leaving a single well pronounced peak 
in the cross section. The standard prescription in the analysis of the unstable top-quark 
production consists of the complex shift E — >■ E + iTf, where E is the top-antitop pair 
energy counted from the threshold 0]. Though this procedure incorporates the dominant 
effect of the finite top-quark width, it does not fully account for nonresonant processes like 
e + e~ —7- tW~b, e + e~ —> bW + t or e + e~ — > W + W~bb where the intermediate top quark is 
not on its (complex) mass shell. Such processes cannot be distinguished from the resonant 
ti production, which has the same final states due to the top-quark instability. Moreover, 
a naive use of the above prescription results in spurious divergences of the cross section in 
the nonrelativistic effective field theory [p8| |30| and it has to be elaborated to make the 
high-order calculations self-consistent. Recently two different approaches have been used 
to refine the analysis of the finite width effect. The first is based on so-called "phase space 
matching" [32] while the second [33] relies on the effective theory of unstable particles 



In particular, in ref. [ 33 1 the NLO nonresonant contribution to the total cross section has 
been computed. 

In the present paper we develop an alternative approach which systematically takes 
into account the effect of top-quark instability. The approach is applicable to the threshold 
production of an unstable particle which is almost degenerate in mass with one of its decay 
products. It introduces into the nonrelativistic effective field theory [^, 36] an additional 
scale p l / 2 rrit characterizing the dynamics of the top-quark decay into a nonrelativistic W- 
boson and an ultrarelativistic bottom quark. The parameter p is related to the difference 
of the top quark and VF-boson masses, p = 1 — mw/™>t, and p 1 / 2 plays the role of the 
heavy quark velocity in standard potential nonrelativistic QCD (pNRQCD) []3^, pq| . The 
new scale obeys the hierarchy vrrit -C p l ' 2 mt < m; or u < p 1 / 2 <C 1 and the cross section 
is constructed as a series in the scale ratios. The method is applicable for the analysis of 
the high order corrections and naturally resolves the problem of the spurious divergences. 

In the next section we outline the main concept and derive the NLO nonresonant 
contribution to the total cross section in the leading order in p. In section |3| we describe 
the calculation of the high order terms of the expansion. The strong coupling corrections 
to the result of section || are computed in section ||. In section || we show how the spurious 
effective theory divergences associated with the top-quark instability are eliminated within 
our approach. Section || collects the final result and numerical estimates. Application of 
our result to the analysis of the experimentally measured cross section with invariant mass 
cuts is discussed in section [?]. Section || is our conclusion. 



2. Finite width effect beyond the complex energy shift 

In the Born approximation the total cross section of top-antitop production in electron- 
positron annihilation is related through the optical theorem to the imaginary part of the 
one-loop forward scattering amplitude shown in figure |l| The corresponding expression 
for the normalized cross section R = a(e + e~ — > tt)/ao, uq = Atto 2 /(3s), in the threshold 
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Figure 1: 



forward scattering diagram corresponding to the leading order top-quark pair 



production process. 



region s ~ Am 2 can be obtained by the standard nonrelativistic expansion of the top-quark 
vertices and propagators in v and reads 

2Q e Q t v e vt , {a 2 e + vl) 



R 
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where the ellipsis stands for the relativistic corrections, Qf is the electric charge of fermion 
/ in units of the positron charge, N c = 3 is the number of colors, and xz = rn 2 z /{Am^) 
with the Z-boson mass m Z - The couplings of fermion / to the Z-boson are 

v f = oT^ > a f 



T 3 



2tS w c w 



(2.2) 



where * is the third component of the fermion's weak isospin and s w (c w ) is the sine 
(cosine) of the weak mixing angle. Note that only the vector coupling of the top quark 
gives the leading order contribution and the axial coupling is suppressed by an additional 
power of v. The last factor in eq. ( |2,lD is 

*^ = /^ra = -^ (2 - 3) 

which is nothing but the Green's function of the free Schrodinger equation at the origin. 
Formally the integral in eq. ( |2.3| ) is linearly divergent but the divergent part is real and 
does not contribute to the cross section. To handle the divergence we use dimensional 
regularization with d = 4 — 2e, where the integral ( |2.3[ ) is finite even for e = 0. The strong 
interaction has a significant impact on the threshold cross section. Close to threshold when 
v ~ a s the Coulomb effects become nonperturbative and have to be resummed to all orders 
in a s by substituting eq. (|2.3j ) with the full Coulomb Green's function 

CFOi s m 2 



G c (0, 0; E) = G (0, 0; E) + G 1 (0, 0; E) 
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where \P is the logarithmic derivative of the Gamma function and Cp = 4/3. The one-gluon 
exchange contribution Gi(0,0; E) is ultraviolet divergent. Again for stable top quarks the 
divergent part is real and does not contribute to eq. ( |2.lD . In the MS subtraction scheme 
this term reads 

CFCH s m 2 
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Let us now consider the top-quark decay. Every decay process is suppressed by the elec- 
troweak coupling constant a ew . We adopt the standard power counting rules a s ~ v, 
OL ew ~ v 2 ■ Thus to NNLO if the top quark decays the antiquark may be treated as a 
stable particle and vice versa. The dominant effect of the top-quark instability is related 
to the imaginary part of its mass operator in diagram |2](a). In the massless bottom quark 
approximation and with the off-shell momentum p the mass operator reads 

lm ^° ){p2)] = T^2 p3 { 1 + 2! f-) (!-|f)V-"40, (2-6) 

where Gf is the Fermi constant and we use the approximation Vtb = 1- Close to the mass 
shell one has p 2 = m 2 — 2(p 2 — mtE) + . . ., where p is the spatial momentum of the top 
quark, and the mass operator can be expanded in z = (p 2 — mtE)/m 2 <C 1 



2 2 



where x = mw /nit, = I* + 0(a s ) and 



+ ... . 

(2.7) 



rf ) = g^(l + 2x 2 )(l-x 2 ) 2 , (2.8) 

is the leading order top-quark electroweak width. Note that in the expansion (|2.7| ) we 
consider 1 — x = p to be of the same order of magnitude as z. The first term in the last line 
of eq. ( |2.7| ) describes the standard shift of the pole position of the top quark propagator into 
the unphysical sheet of the complex energy plane characteristic for unstable particles. After 
Dyson resummation it replaces the argument of eq. fl2.3| ) by E + iTt , which is the original 
prescription of ref. [Q]. Eq. ( p.7| ), however, has the remainder which also contributes to 
the imaginary part of the forward scattering amplitude. Since the remainder vanishes for 
on-shell top quark, it represents the nonresonant process e + e~ — > bW + i or e + e~ — > tW~b. 

In the nonresonant contribution the integral over the virtual momentum p is saturated 
by the region |p| ~ p l l 2 mt. The main idea of our approach is that if p is considered as 
a small parameter this momentum region corresponds to a nonrelativistic top quark with 
the energy pQ — mt ~ p 2 /rrit ~ pmt, and one may apply the well elaborated pNRQCD tools 
for the calculation of the cross section. In this complementary nonrelativistic expansion 
the heavy quark velocity v is replaced by p 1 / 2 and we have the hard scale mt, the soft 
scale p 1//2 mt, and the ultrasoft scale pnif. Note that in this case the VK-boson is also 
nonrelativistic while the bottom quark is ultrarelativistic with the ultrasoft four- momentum 
of order pmt 3> m^. Further expansion of eq. (^.7|) in p ~ z gives 



ImP (0) (*)] = y-y 



e(z-p) + [j-^)0(p-z) + O(p,z) 



(2-9) 
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Figure 2: 



forward scattering diagrams containing bW + t and tW b cuts. 



By inserting the second term of eq. ( |2.9| ) into the diagram g(a) one obtains the following 
contribution to the nonresonant cross section 
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where 5\ 



The leading order result reads 
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In eq. ( gH| ) we neglect the ra(£ term in z since it gives a subleading contribution sup- 
pressed by mtE/p 2 ~ v 2 / p <C 1 according to our scale hierarchy. Note that one can obtain 
the result ( 2.1 1|) without subtracting the first term of eq. fl2.9|) by direct evaluation of the 
corresponding Feynman integral in dimensional regularization. In this case the resonance 
contribution vanishes since eq. (^^) becomes a scaleless integral. 

Eq. ( 2,1C| ) however does not give the full nonresonant contribution and one has to take 
into account all diagrams of this order with the bW + i or tW~b cuts given in figure The 
calculation is significantly simplified within the nonrelativistic effective theory where all the 
propagators which are off-shell by the amount mt collapse, giving rise to new effective theory 
vertices. From the practical point of view, however, it is more convenient to directly expand 
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the full theory Feynman integrals in p than to use the effective theory Feynman rules. We 
found that beside diagram 0(a) only diagram |2](g) gives a leading order contribution in p. 
After the expansion the off-shell neutrino propagators shrink to points and the diagram 
becomes similar to figure ||(a) . In this case the contribution comes from the imaginary part 
of the vacuum polarization operator of the heavy-light vector current correlator 
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where p 2 
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2(p 2 — mtE) + . . .. It has the following expansion 
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The corresponding contribution to the nonresonant cross section reads 
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Each VF-boson propagator brings a factor (2pm 2 — p 2 ) 1 to eq. ( 2.15 ) which is regular at 



p = so that the diagram does not have a resonance contribution. Finally for the leading 
nonresonant contribution to the cross section we get 



Rr, 




2 2 2Q e Q t v e v t (a* + v 2 e )v 2 



1 



9\/2 
~32" 



In (I + V2 




(2.16) 



Let us now compare our approach to the one of [33]. In ref. |33| the scales mt and p l / 2 mt 
are considered to be of the same order and are integrated out simultaneously. As a result 
the nonresonant contribution is represented by the imaginary part of the Wilson coefficient 
of the local four-fermion e + e~e + e~ operators, i.e. every diagram in figure [2] shrinks to a 
point. The method requires more complex calculations but gives the exact dependence of 
the cross section on p. However, in the next section we show how to compute a sufficient 
number of terms of the expansion in p to ensure very good accuracy of the approximation 
for the physical value p ~ 0.53. At the same time application of the method |^|] to the 
calculation of the strong interaction corrections to the nonresonant contribution seems to 
become technically complicated while our approach does not as we show in section ^. 
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3. Relativistic corrections 



To obtain the high order terms of the nonrelativistic expansion of the nonresonant cross 
section in p we use the method of regions [39, [EJ. By the optical theorem the problem 



is reduced to the calculation of the contribution from bW + i and tW~b cuts to the imagi- 
nary part of the two-loop electron-positron forward scattering amplitude, figure As an 
example let us describe the evaluation of diagram §(e). The corresponding scalar integral 
reads 

d d k d d l 

(I 2 + ie)(k 2 + q-k + ie)(k 2 - q ■ k + ie)[(k + I) 2 + q • (k + I) + m 2 p(2 - p) + ie] 

' (3-1) 



[{k + l) 2 -q-{k + l) + m 2 p(2 - p) + is] ' 

where q = (2mt,0) is the photon/Z-boson momentum corresponding to the top-antitop 
threshold. In the limit p — ¥ the only nonvanishing contribution to the imaginary part 
comes from the region of potential momentum k and ultrasoft momentum I 

ko ~ rritp , k ~ rritp 1 / 2 , I ~ uitp . (3-2) 

By imposing this scaling we expand eq. (3.1) in p. For example in the leading order we 



obtain 



(I 2 + ie)(-k 2 + 2m t k +ie){-k 2 - 2m t k + ie)[-k 2 + 2m t (k + l ) + 2m 2 t p + ie] 
2—^- ( 3 - 3 ) 



[-k z - 2m t (k + Iq) + 2m\p + ie] 

Then we evaluate the integrals over the zero components of the loop momenta by closing 
the integration contours in the upper half of the corresponding complex planes. At this 
step we have to distinguish between the different cuts of the diagram and pick up only those 
poles which correspond to bW + i and tW~b cuts. The integral over k becomes a one-loop 
massive tadpole integral in d — 1 dimensions and can be performed easily. Integration over 
I yields a Gauss hyper geometric function 2F1 with half-integer parameters, which can be 
expanded in e with the Mathematica package HypExp 2 |jH|| . 

The apparently more complicated diagrams |](g)-|2]C/) with t-channel neutrino propa- 
gators do not pose any new problems since the off-shell neutrino propagators do not depend 
on the loop momentum after the expansion and the resulting integrals can be computed 
in the same way. We use QGRAF [^], q2e, and exp [[E|, to generate the diagrams and 
produce FORM-readable expressions for the amplitudes. The expansions and integrations 
are performed with custom code written in FORM [J45|| . The result of the calculation is 
presented in section |6[ 

4. Strong coupling corrections 

In the leading order in p the a s corrections are obtained by the gluon dressing of dia- 
grams |2](a) and 0(g) shown in figure ||. Since the top quark is nonrelativistic we can apply 
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Figure 3: e + e forward scattering diagrams contributing to the 0(a s ) correction in the leading 
order in p. The curly lines denote gluons. 



the standard pNRQCD arguments to identify the relevant regions of virtual momentum. 
For diagram |3](a) the leading contribution comes from the Coulomb gluon with the po- 
tential momentum qo ~ q 2 /mt ~ pmt- As in the case of the Coulomb gluon exchange in 
pNRQCD, the corresponding correction is power enhanced due to the Coulomb singularity 
but the enhancement factor here is 1/p 1 ^ 2 rather than 1/v. An important difference with 
respect to pNRQCD is that since a s j p 1 ! 2 <C 1 one does not need to resum the Coulomb 
corrections to all orders. Following the analysis of section |2] one gets the potential Coulomb 
gluon contribution of the following form 
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(4.1) 



where = In f \ an d 6^ defines the correction to eq. ( 2.1Q|) , 5\ = 5^ + 5^ + 

0{a 2 ). Note that one has to keep a nonzero (complex) energy in the top-quark propagator 
in eq. ( |4.1| ) since it serves as an infrared regulator. The contribution of the hard gluon 
with momentum q ~ mt in diagram |3](a) can be related to the Wilson coefficient in the 
nonrelativistic expansion of the vector current j M 

d v 



j = Ou^trx + -A I ip Jf crD 2 x + 



(4.2) 



where tp and x are the nonrelativistic quark and antiquark two-component Pauli spinors 

and the Wilson (matching) coefficients are c v = 1 — 2CFa s /ir + . . . and d v = 1 + The 

0(a s ) term in c v results in the following correction to 5\ in eq. (|2.3D 
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Diagram ||(6) is determined by the correction to the mass operator in the limit z — > that 
can be read off the result of ref. 
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The correction factor should in principle be substituted into the integral ( 2.11 ), However 
all the terms in eq. fl4,4|) except the last logarithmic one actually describe the corrections 



to the top-quark width 



_|_ pj 1 ^ w here 
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They all disappear when the leading order top-quark width in eq. (2.10) is replaced with 
its corrected (physical) value. The remaining term gives 
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The same procedure applies to diagram ||(<i) which is determined by the correction to the 
polarization function 
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After factoring out the corrections to the top-quark width one gets 
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where Li 2 stands for the dilogarithm function. 

The potential gluon contribution to the nonfactorizable diagram ||(c) in the leading 
order in v vanishes for the total cross section |29]. The hard gluon contribution to this 
diagram is power suppressed. Thus the diagram vanishes in our approximation. This 
completes the calculation of the dominant 0(a s ) corrections. In the next section, however, 
we present the analysis of the corrections suppressed by an additional power of p 1 / 2 which 
addresses an important issue of the spurious divergences in the nonrelativistic effective 
theory of unstable top-quark production. 



5. Eliminating the spurious divergences 

Let us first outline the problem. In the pNRQCD perturbation theory the Coulomb Green's 



function gets corrections due to the second term of eq. (^2) and from the relativistic 
corrections to the Coulomb Hamiltonian 
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of the following form 
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3 m t 



(5.2) 



As it has been pointed out the Green's function in the above equation is divergent and in 
dimensional regularization reads 
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where the standard MS factor 



An 



per loop is suppressed. For real energy values the 
divergent part of eq. (5.2) is real and does not contribute to the cross section. After the 
complex energy shift the divergence gets the imaginary part proportional to I\ 
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resulting in a divergent cross section. In the previous analysis the Green's function was 
renormalized, as in eq. ( |2.5| ), leaving a finite but scheme-dependent cross section. The 
solution of this problem is straightforward within our approach. Indeed the logarithmically 
divergent integral in eq. ( |5.4[ ) has a physical cutoff scale /5 1//2 mj where the imaginary part of 
the top-quark mass operator vanishes. In the above expression this scale is set to infinity. 
Thus within the expansion by regions framework eq. ( |5.4| ) represents a contribution of 
the pNRQCD potential momentum region \p\ ~ \fm%E. To get the total result one has 
to add the contribution of the additional potential region \p\ ~ p l l 2 mt. The latter is 
given by the interference of the corrections to the Green's function due to the top quark 
mass operator (E 



and due to the relativistic corrections to the vector current (4.2) and 
the Hamiltonian (|5.1[ ) in the second order of time-independent perturbation theory. The 
resulting integral after expansion in E/(pmt) is infrared divergent and in dimensional 
regularization reads 
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In the sum of eqs. (JlO|) and (^]5|) the poles in e cancel each other and one gets the finite 
result for the correction 
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One formally reproduces the sum of the two regions, eq. ( |5.6[) , if in eq. (5.2) the MS 
renormalized Green's function ( |2.5[) is evaluated with p = e~ 5 ^p 1 ^ 2 mt, which can be taken 
as a practical prescription for the calculation of the cross section. A similar analysis in the 
case of P-wave heavy-quarkonium production has been performed in refs. [28, pCj] . Note 
that eq. ( |5.6D gives only a part of the 0(p l / 2 a s ) corrections to the nonresonant cross section 
which corresponds to specific terms in the nonrelativistic expansion of diagram ||(a) and 
does not account for the axial coupling of the top quark and the higher order terms in the 
expansion of the mass operator as well as the contribution of other diagrams which vanish 
in the lower orders. The complete result for the 0(p 1 ' 2 a s ) corrections can in principle be 
obtained within the approach described in section [|. However, as we will see in the next 
section, the leading term of the expansion in p gives a good approximation of the total 
series and is sufficient for the practical applications. 



6. Results 

The NLO nonresonant contribution to the threshold top-antitop production in electron- 
positron annihilation reads 



8iV c Ft 
Trp 1 / 2 m t 

2 2 2Q e Q t v e v t {a\ + v 2 )v 2 



l-xz 



cVV 



+ 



(«e + V 2 e )a 2 t AA 



+ 



+ 



+ 



+ 



2 2Q e v e c u 



Q'i + 



QlQl + 



(l-x z ) 2 
2Q e Q b v e (a b + v b ) 



2Q e Q t v e a t 2(a 2 e + v 2 )a t v t 



l-x z 



+ 



l-x z ) 2 
a'i + v 2 )(a b + v b ) 2 ' 

Jb 



l-xz 

,2 i „,2 



J a 



VA 



+ 



{a 2 + v 2 e )c\ 



QeQtQb + QeV, 



X Z ) s^(l-x z 
Qt{a b + v b ) + Q b v t 



fc + 



Q e Q b v e a t 



+ 



\-xzY 
(a 2 + v 2 )a t (a b + v b ) 



l-xz 
(flfe + v b )(a 2 e + vl)v t 



l-xz 



QeQt + QeV e 7T. T 

S W {1 - X Z ) Hi 

rVn, , ^ QbCw + (a b + v b )s 



w (a 2 + v 2 e )v t c u 



,(l-x z ) 



(l-xz) 



xz 



(a 2 + v 2 )a t c u 



l-xz s w (l-x z )' 



+ 



s w (l - x z 
QeQt (a e + v e )vt 

s w s w(l ~ x z) 

QeQb , (a e + v e )(a b + v b ) 



w (aj + v 2 )(a b + v b )c v 



+—f 9 + 



+ 



4(1 - x z ) 



s w (l - x z f 

y (gg + Ve)at , A 

-x z ) 
(a e + v e )c u 



ff 



st,{l 



fi + 



9± 



[l-xz) 



fj 



(6.1) 



where /„ = f n (p) stands for the contribution of diagram ||(n) and the superscript distin- 
guishes the contribution of the vector and axial top-quark coupling. Each function f n is 
found as a power series in p up to 0(p 12 ). A text file with the expressions for all f n is 
attached to the paper's source files on the arXiv. The prefactor of eq. (|6.1| ) is chosen in 
such a way that fY v = l + 0(p). In order to cross-check our results we also expanded the 



integral representations given in rcf. [33) and found perfect agreement 
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Figure 4: (a) Dashed lines represent the expansion of the function fj(p) in p through 0{p ) for 
n = 1, . . . , 6. (6) Dashed lines represent the [1/11], [11/1], [9/3], [6/6], and [3/9] Pade approximants 
of the series for the function f¥(p). In both pictures the vertical line marks the physical value of 
p and the solid curves represent the all-order numerical result of ref. |33| . 



With the exception of and f£, the series converge very well for the physical value 
p ~ 0.53. As an example, in figure ||(a) we compare the expansion of the function fj with 



the exact numerical result obtained from the integral representation given in ref. [33]. For 
the two special cases we perform a Pade resummation to improve the convergence, i.e. we 
construct the Pade approximants 



[n/m] 



En i 
i=k a iP 

i + TZi M 



(6.2) 



where k equals 1 and 2 for f% and f^, respectively. The coefficients Oj and b{ are determined 
by matching the expansion of eq. fl6.2|) in p to the series for and fj~. In figure |||(&) 



we compare different Pade approximants of f¥ with the result of ref. |33| and find perfect 
numerical agreement, which means that Pade resummation solves the problem of slow 
convergence. For our numerical analysis of the cross section we choose the approximant 
[9/3] for f¥ and [8/3] for f£, which deviate from the result of ref. |£| by less than one 
percent for the physical value of p. 

The NNLO nonresonant contribution in the leading order in p is given by the 0(a s ) 
correction to eq. (|2~Trj ). By combining eqs. flO]) , flO|) , (gig ), and (|J) one gets 
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(6.3) 
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Figure 5: The NLO nonresonant contribution to the cross section (p.l|) (solid curve) and the 



leading order approximation ( 2.16 ) (dashed line). The vertical line marks the physical value of p. 



Let us now study the numerical effect of the correction. We adopt the following input 



values []48| 



m t = 172 GeV , m w = 80.399 GeV , m z = 91.1876 GeV . 
G F = 1.16637 • 1(T 5 GeV" 2 , a s (m z ) = 0.1184 . 



(6.4) 



In eq. ( |6.3D we use the value of the strong coupling constant a s (/J,) = 0.1129 corresponding 
to the physical normalization scale [i = p l l 2 mt. This value is obtained from a s (mz) by 
means of the RunDec program [fl9[|. 

Our total result for the NLO contribution ( |6.1D is plotted in figure [| as function of p 
along with the leading order term (2.16). The latter turns out to be a good approximation 
in the whole interval < p < 0.6 and deviates from the total result by less than 5% 
at the physical value of p. This justifies our approximation of the NNLO nonresonant 
contribution fl6.3|) by the leading order of the expansion in p. The numerical effect of 
the nonresonant contribution on the total threshold cross section is shown in figure ^. In 
this plot we use the leading order pNRQCD approximation for the resonance contribution 
corresponding to the Coulomb Green's function (|2.4j ) with the strong coupling constant 
normalized at the soft scale pL s = a s (p, s )CFmt. 



7. Measured cross section and invariant mass cuts 

Experimentally the top-antitop pairs produced in electron-positron annihilation are to 
be reconstructed from the lepton plus four jets final state or all hadronic six jet events. 
Realistic simulations show that the invariant mass p 2 = {jpw+Pb) 2 OI the three jets resulting 
from the bottom quark and the hadronic decay of the VK-boson can be determined and 
the corresponding invariant mass distribution can be obtained with high accuracy ||50| . 
The total cross section of the threshold production considered in this paper is obtained by 
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Figure 6: The normalized cross section of top-antitop production in electron-positron annihila- 
tion as function of the energy counted from the threshold. The dotted curve represents the leading 
order pNRQCD Coulomb approximation. The dashed curve includes the NLO nonresonant contri- 
bution and the solid curve includes the NNLO nonresonant correction as well. No strong coupling 
corrections to the resonance contribution are included. 



integrating this distribution over the kinematically allowed interval m 2 ^ < p 2 < mf , which 
corresponds to the integral over p in eqs. ( 2.11 ) and ( |2.15 ). One may suggest that a tight 
invariant mass cut m 2 — p 2 < A 2 with A ~ (ratT^) 1 / 2 would separate the production of 
"true" on-shell top quark and antiquark from the nonresonant background. However, it 
is not possible to suppress the nonresonant contribution without significant modification 
of the resonant one. The resonant contribution of unstable top quark corresponds to the 
Breit-Wigner shape of the invariant mass distribution, which falls off rather slowly as the 
invariant mass deviates from rut and is strongly affected even by loose cuts. To quantify 
this statement let us consider a loose cut (mjTt) 1 / 2 <C A. In the Born approximation it can 
be implemented by replacing the physical cutoff pm 2 in the argument of the theta-functions 
in eqs. fl2.11| ) and fl2.15| ) with A 2 /2. For A <C p x l 2 mt one gets 

*f--T + °fTr-) (7J) 

7r A \p 1 ' A m t J 

and #2 ~ (A/ / o 1 / 2 m t ) 3 . The leading term in eq. (7.1) corresponds to the first term 
of eq. ( |2.11|) and represents the modification of the resonant contribution to the cross 
section ( |2.1| ) by the cut. It is merely the dominant effect of the cut unless A ~ / o 1//2 mj. 
The dependence of the correction to the cross section on the invariant mass cut value is 
shown in figure ^. Starting from A ~ 80 GeV it is relatively weak and the correction is 
well approximated by eq. ( 2.16| ). This gives a low bound on the acceptable cut value which 
can be used for the determination of the total cross section. One may study the cross 
section with tighter cuts as well. However all the high order QCD results for the total 
cross section are not applicable in this case and accuracy of the theoretical predictions 
would be significantly reduced. 
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Figure 7: The nonresonant contribution as function of the cut on the three-jet final state invariant 



mass. The horizontal line corresponds to eq. ( 2.16| ). 



8. Summary 



We have developed a new method for the analysis of the top-quark instability in the 
threshold top quark-antiquark pair production in electron-positron annihilation based on 
the nonrelativistic expansion in the parameter p = 1 — mw/ m t- Within this framework we 
obtain the NLO nonresonant contribution to the total threshold cross section overlooked 
in the standard analysis and confirm the result of ref. [33] obtained within a different 
effective theory approach without the expansion in p. The NLO contribution is negative 
and amounts to about 3.1% of the leading cross section above the threshold and competes 
with the LO contribution below the resonance region. We extend the analysis to the 
NNLO 0{a s ) nonresonant contribution which is computed to the leading order in p. The 
corrections involve In p terms which are a new type of the large logarithms in the theory of 
top quark threshold production. The NNLO nonresonant contribution amounts to —0.9% 
at the threshold. 

Our method can also be applied to the calculation of the nonresonant contribution to 
the threshold top quark-antiquark pair production in hadronic collisions. The process is 
of particular interest since a significant number of top quark-antiquark pairs is going to be 
produced at the LHC near the threshold and the accuracy of the top quark reconstruction is 
expected to be sufficiently good to study the threshold region. A comprehensive analysis of 



the process is given in refs. p5\ , 52], where the NLO approximation without the nonresonant 
part has been used to derive the top-antitop invariant mass distribution. The missing 
nonresonant contribution can be directly obtained from our result. Indeed, in the threshold 
region the cross section is dominated by the color singlet top-antitop configuration which 
is produced mainly through gluon fusion. Thus, the dominant nonresonant contribution is 
entirely due to the imaginary part of the top-quark mass operator and in full analogy with 
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eqs. (|27T[) and (|2.10|) one gets 



^ nr (hadrons ^ tt) = — f — - ) i?? e T n (hadrons ->• tt) L_. r (1 + £>(p, a s )) , 
37r \pm t J II-. -/j ; - 



which should be added to the resonant contribution considered in refs. pi 
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